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ABSTRACT 


The purpose of the Thrust Vector Control (TVC) System is 
to point the engine thrust at the vehicle center of mass and thus 
maintain attitude stability during propulsive maneuvers. This is 
accomplished by mounting the engine in a two-axis gimbal system. 
The TVC system then controls the pointing of the engine by closed 
loop control of two linear actuators which extend or retract and 
rotate the engine in its gimbal system. 

The velocity vector (AV) pointing error incurred during a 
propulsive maneuver is made up of several parts, and the TVC 
system error is one of these. This error arises because the 
center of mass and engine thrust are not perfectly aligned. 

It is the purpose of this memorandum to assess the magni- 
tude of the TVC system portion of the AV pointing error and to 
provide simplified and linearized models for predicting its magni- 
tude for the various propulsive maneuvers. The accuracy of the 
linearized analysis is verified using a detailed nonlinear six- 
degree -of-freedom computer simulation. 

The Jet Propulsion Laboratory is responsible for the 
Viking Orbiter System, which is part of the overall Viking Project 
managed by the Viking Project Office at Langley Research Center 
for NASA. The spacecraft will be launched on a Titan III-E/ 
Centaur launch vehicle in August 1975. 
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I. INTRODUCTION 


The purpose of the thrust vector control (TVC) system is to point the 
engine thrust at the spacecraft center of mass (CM!) and. thus maintain attitude 
stability during propulsive maneuvers. Due to various errors, the engine 
thrust vector will, in general, be moved from its initial position, and the 
spacecraft velocity change (AV) will not be exactly in the direction com- 
manded. The total AV error is made up of several components, e. g. , com- 
manded turn errors, gyro drifts, etc., and the TVC system error. The TVC 
system error arises because CM and the engine thrust are not perfectly aligned. 
This alignment error is comprised of two fundamentally different error sources. 
(1) CG error, which arises because the position of the actual spacecraft CM 
and the position of the actual engine center of thrust are not known with 
complete accuracy; and (2) angular error, which comes from several sources 
and results from the fact that the direction of the actual engine thrust vector 
is not known with complete accuracy. 

It is not the purpose of this memorandum to assess the magnitude of 
these alignment errors, but instead to predict the TVC system contribution 
to the total AV pointing error as a function of these alignment errors. In 
other words, given so much CG error and so much angular error, the analy- 
sis presented here will enable the prediction of the TVC system AV pointing 
error. 

The analysis presented herein is based on the work originally done dur- 
ing the design of the Mariner Mars 1971 thrust vector control system. The 
differences are principally due to spacecraft structural dynamics and mass 
properties, propellant slosh dynamics, and the TVC system hardware 
implementation. 
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II. POINTING ERROR DEFINITION 


As previously mentioned, the AV pointing error is comprised of such 
errors as gyro-drift error, commanded turn error, and TVC system error. 
The TVC system error consists of transient and steady-state characteristics. 
Before developing the TVC system pointing error characteristics, a brief 
discussion of the basic system operation is necessary. A sketch of the 
Viking Orbiter 1975 (VO 1 75) engine gimbal system is shown in Fig. 1. The 
thrust pointing error can be expressed as 

v = 6 + "cm (1) 


whe re 

6 = inertial position of the spacecraft as sensed by the inertial 

reference unit (IRU) 

&CM = entire angular position uncertainty of the engine CM system 
The actual physical motion of the gimbal can be expressed as 


Q 


G 


a 


CM 


- a 


( 2 ) 


where 4 total angular error of the engine - gimbal system and including 

any pre-aim null shift. 

Denoting the TVC system steady- state gain by K, 

a G = K8 


From Eqs. (1) and (2) 


<P 




CM 



+ CX 


CM 


( 3 ) 
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Inspection of Eq. (3) shows that a TVC system with a high DC gain will 
exhibit a steady-state error that is insensitive to angular errors but allows 
CM shift (positional) errors to pass unattenuated. In the Viking Orbiter 
(VO' 75) vehicle, as in the case of Mariner Mars 1971 (MM’ 71) vehicle it is 
known that the CM will undergo large positional changes so that i- s ^h.e 

dominant source of error. For this reason, a ’’trim loop" or "path guidance 
loop is employed which, in effect, sets K equal to - 1 and Eq. (3) becomes 


III. TRANSIENT RESPONSE 

A. Dynamic Model 

To obtain an accurate estimate of TVC system performance, especially 
during the initial transient, a complete nonlinear simulation must be per- 
formed. This has been done, but it is unduly expensive for general accuracy 
determination. Instead, the TVC system accuracy is determined using linear- 
ized models, and the results are checked using the detailed simulation. 

Fig. 2 shows a simplified block diagram of the TVC system. Here the gyro 
and integrator in the IRU are modeled by gains only and then mixed with a 
rate -to-position time constant of T^,, The forward compensator is a gain 
and path guidance is a single positive feedback lag. All the other lags 
in the system are lumped into a single time constant T^ associated with the 
gimbal servo. The torque on the spacecraft is -K^-, (4> - q q) which acts upon 
a simple double integrator approximation (rigid-body) of the attitude dynamics. 

B. Velocity Vector Error 

Equation (1) gives the steady- state thrust pointing error. Using Eq. (2), 
the following form is obtained: 

<p - 6 + a + a (5) 

This equation is true, in general, and otq is the physical motion of the gimbal 
system. Thus, if <|>(t) is defined as the time variation of gimbal angle, cp 
can be expressed as 
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<p(t) = e{t) + <|>(t) + <* A ( t) 


( 6 ) 


as the general equation. 

Let the mass of the vehicle be M, and the magnitude of the engine 
thrust be T, and assume that these are constant over the duration of the 
transient. Further, assume small angle nature for 4>(t) and neglect second- 
order terms. Under these assumptions, the magnitude of the velocity change 
is 


AV(t) 



(7) 


and the component of AV normal to the desired direction (side velocity) is 
given by 


AV x (t) 


r_ 

M 


(p{r) dT 


( 8 ) 


The pointing error of the AV vector is thus approximated by 


Y(t) 


AV x (t) 
A V(t) 


(9) 


If AV^(t) is known, y(t) can be determined by small angle approximation. 

Denoting the independent variable of the Laplace transformation by s, 
Eq. (6) becomes 


<P( s) = 8(s) + 4>(s ) + a A (s) 


and Eq. (8) is 


AV x (s) 



( 10 ) 


(H) 
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Defining 


H(s) 


A 6( 5 ) 

“g ( s) 


and using Eq. (2) in transform form yields 

<p(s) = H(s ) Q' cm ( s ) + <M S ) + [l “ H (s)] <>' A ( S ) 


( 12 ) 


( 13 ) 


Further define 


G(s) 


A ^(s) 

a G ( s ) 


(14) 


then 


tp(s) = [H(s) + G(s)] or CM (s) + [l - G(s) - H(s)] ff A (s) (15) 


and Eq, (11) becomes 


AV x (s) 


_ X fH(s) + G(s)~ l „ , , + [I 

' M [ s J CM (S ' + L 


- H(s) - G(s ) 
s 


K 


(s) 


( 16 ) 


As expected with a linear system, superposition holds and CM positional off- 
set errors can be treated separately from the angular errors. 

C. Transfer Functions 

Referring to Fig. 2, the previously defined transfer functions can be 
determined : 


H(s) 


= 6(s) 

' “g ( s) 
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After some manipulation. 


whe r e 


and 


K 


FG 


H(s) = 


T A T P 2 T P + T A (1 " K PG ) 

s H : s + 1 


1 - K 


PG 


1 - K 


PG 


D(s) 


(IV) 


D(s) = 


T T J 
P A 


k t k g k fc k a 


4 , 
s + 


J Fa - K PG> + T 


K T k g k fc k a 


J (1 - K ) 

rp rn i 

P G K 


T “G “FC K A 


s + (T g + T p ) s + 1 


(18) 


1 - K 


K 


PG 


FG 


k g k fc k a 


(19) 


K 


PG 


K 


FC 


K. 


( 20 ) 


G(s) was defined as 


G(s) 


4 > ( s ) 
o G ( s ) 


and from Fig a 2 


G(s) 


(T p s + 1) (T g s + 1) 

D(i) 


( 21 ) 


D» Numerical Values 

Table 1 gives the values for the parameters for the transient response 
model for the first midcourse maneuver, first orbit trim, and second orbit 
trim (post lander release) conditions. As previously mentioned, the TVC 
system DC gain must be - 1 in value; thus 
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( 22 ) 


k pg " 1 4 k g k fc k a 

Equation (22) is not an exact equality in Table 1, because nearest standard 
values are used for the TVC select resistor (Ref. 3) which set the forward 
compensator and path guidance amplifier gains. The last entries in Table 1 
are for the torque constant Krj, given by the engine thrust-lever arm product. 
When the first midcourse values are substituted into Eq. (16), side velocity 
due to CM offset error is 


AV X^ S * _ 10. 406 (1.217 s + 1) 

0. 5762 s 4 + 3. 446 s 3 + 7. 086 s + 6. 860 s + 1 


(23) 


and for angular error 


AV X^ _ 1.274 (0. 5762 s 4 + 3.446 s 3 - 2. 853 s 2 - 1. 308 s + 1) 
s (0. 5762 s 4 + 3.446 s 3 + 7. 086 s + 6. 860 s + 1) 

In a similar manner, when the first orbit trim values are substituted, 


AV X (S * _ 15. 331 (1.217 s + 1) 

0. 3537 s 4 + 2. 116 s 3 t 8. 103 s 2 + 6. 860 s + 1 


(25) 


and for angular error 


AV X^ _ 1. 877 (0. 353 7 s 4 +'2, 1 16 s 3 - 1. 836 s 2 - 1. 308 s + 1) ^6) 

s (0. 3537 s 4 + 2. 116 s 3 + 8. 103 s 2 + 6, 860 s + 1) 


and finally for the second orbit trim condition, 


AV X* s) 41. 391 (0. 9345 s + 1) 

a CM^ 0. 5852 s 4 + 3. 768 s 3 + 8. 833 s 2 + 6. 860 s + 1 
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and for angular error 


AV x (s> 

°’A (S) 

These equations are mathematically more tractable when put through a partial 
fraction expansion. For each condition, a partial fraction expansion will yield 
equations of the following form: 


3, 726 (0, 5852 s 4 -f 3, 768 s 3 - 1, 547 s 2 - 4, 248 s + 1) 
s (0, 585Z s 4 + 3. 768 + 8, 833 s 2 + 6, 860 s + 1) 


(28) 


and 


AV x (s) 

a CM^ 


K 


K. 


s + a s + b 


Z 

s + w 


Z 

s + w 


(29) 


AV x (s) 

°A (S) 


K. 


K. 


Ms 


s+a s+b s+w s+w 


(30) 


whe re 

K^, K^, a, b = real constants 

Z,w = complex constants 

Z, w = complex conjugate of Z,w 

Taking the inverse Laplace transform of Eqs. (29) and (3 0) will allow the 
calculation of AV" x (t) for any CM offset and angular error for each -of the 
flight conditions. The AV pointing error can then be calculated using Eq, (9), 
In the time domain, Eqs, (29) and (30) become for unit CM and angular errors 
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( 32 ) 


K 1 / -at \ K 2 /, -bt \ 

AV XA (t> = AV(t) l 1 ' 5 )-~V ^ ) 

^(l-e' Wt ) + | (I-''™) 

The constants in Eqs, (31) and (32) are listed in Table 2 for each of the 
three flight conditions. 

IV. STEADY -STATE ACCURACY 

In developing the relationships for transient response error, all param- 
eter values were taken as nominal. This is, of course, the ideal case and 
with the actual hardware there will be not only tolerances on the gains but 
cross -coupling of the axes. It is the purpose of this section to estimate these 
factors and predict the steady-state AV pointing error due -to the TVC system. 

A. Position Gain Equations 

Figure 3 shows the position signal path for one axis (pitch or yaw) of 
the TVC system. Denoting X as pitch and Y as yaw, the pitch and yaw com- 
pensator steady-state outputs are 


K GY k fcx 


AX 


1 - K 


PGX 


’X 


(33) 


V 


k gy k fcy 


AY 


1 - K 


(34) 


PGY 


From Fig. 3, it can be seen that the extensions of the gimbal actuators are 


e 


X 


V AX K PV + V PAX K PA 

k fbx k fbx 


(35) 


V AY K PV V PAY K PA 

k fby k fby 


(36) 
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where Vp^ and Vp^y are the pre-aim voltages, and if this bias is 
separated out, 


e x 6 X + e PX 


^ y "t e p y 


whe re 


Y K 
AX DV 

K FBY 


V K 
AY DV 

K 

FBY 


and the pre-aim extensions 


V K 
PAX PA 

k fbx 


V K 
PAY PA 

k fby 


Combining Eqs, (33), (34), and (39), (40) 


■ M 


with 


6 — col "I 6y | 


e = col p x , e y } 
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and 


K FD 





(44) 


where 



k gy k fcx k dv 

(i - Kp GX ) K fbx 


(45) 



k gy k fcy k dv 

(! - K PGY^ K FBY 


(46) 


B. Selection of Position Gains 

As previously described, the gains are selected for proper path guid- 
ance, producing steady- state positions 

6 X = " °GX 


e 


Y 


a 


GY 


G x and 0 y are the attitude rotations about the pitch and yaw axes and 

o and a are the components of actual gimbal rotation. For the purpose 
GX GY _ 

of gain selection, the gimbal system is considered linear and that 



(47) 

(48) 
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with being the lever arm from the actuator clevis to engine center of 
rotation. Combining Eq. (43) with the approximations (4 7) and (48) and further 
assuming that the gimbal system is exactly aligned with the spacecraft XY Z 
system, yields 



(49) 


To satisfy the path guidance criterion, it is necessary to select 

K X = K Y " l e (so) 

in order to obtain 


— — 


— — 

O' 


0 

G x 


X 

OC r 

g y 


e Y 

— 





C. Gain Accuracy 

The basic TV C loop mechanization is shown in Fig. 4, Because of the 
large change in vehicle moments of inertia that occur at lander separation, 
compensator and path guidance gain changes (via CCS command) are required 
to maintain system stability. The TVC DC gains as well as the transient 
response are set using "select resistors" that are located for easy access on 
top of the modules and on the multilayer board. These resistor values are 
computed and rounded to nearest standard value by a special computer pro- 
gram. The TVC select resistor values for both flight spacecraft and the proof- 
test model (PTM)/ spare are the same and are listed in Table 3. Using the 
subscript H to denote high gain (pre-lander release) and L to denote low 
gain (post lander release), the compensator and path guidance loop gains are 
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(52) 


K 


FC. 


[ R 12 (R 19 + R Z3 ) + R 1 9 R Z3J 

R 0 -2 R/ 

2 3 6 


R 


K 


12 t +R 19( R 23 + R 4 + r1 1z). 


FC 


H 


R 


(53) 


K 


R 6 R 14 K FC l 
PG L = < R 16 + R 17 )R 9 


(54) 


R 14 R 6 K FC 


K 


H 


PG 


H 


R 


R 16 R 9 2+ R 


16 


(55) 


For assessing the accuracy of the gains, the expedient will be used that all 
error sources are gaussian and independent and worst-case values will be 
taken as 3<r. The worst-case tolerance on resistor values is taken as ±2 per- 
cent, This is made up of an initial tolerance of ±1 percent, a variation due to 
temperature (100 ppm at 50° C) of ±0. 5 percent and a life change of ±0. 5 per- 
cent, The gain variation from a ±2 percent resistor tolerance was computed 
for each resistor in Eqs. (52) through (55) and the results root-sum-squared 
to obtain the total 3cr variation. The results are summarized in Table 4 and 
include the quoted variations in the IRU position scale factor and the 

gimbal actuator feedback scale factor K^. The largest worst-case varia- 
tion is in K_, 0 and is due to the "soft" supply which generates the LVDT 
r io 

excitation voltage and its variation with LVDT input impedance, 

D, Gimbal Transfer Characteristic 

The actual gimbal system transfer characteristic which relates gimbal 
angle and actuator extension is nonlinear and includes inter -axis coupling. 
This transfer characteristic was originally derived for the MM'71 gimbal 
system and was repeated for the VO' 75 system and the results were used for 
developing a computer program for determining the computer command 
subsystem (CCS) pre-aim commands and for developing second-order 
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approximations for the accuracy analysis- To determine the required gimbal 
pre-aim, it is necessary to express actuator extensions as a function of gim- 
bal angles. This is done by expressing first the given gimbal rotations in 
terms of Euler angles and then considering the geometrical constraints in the 
rotated system- The result is a nonlinear system of the form 


X - g x ^X’ 

4>y) 

(56) 

Y = g Y ^X* 

+Y> 

(57) 


However, the reverse problem, to find the gimbal angles when the actuator 
lengths are given, can only be solved numerically. For the purpose of the 
accuracy analysis, approximations can be derived. Referring back to 
Eqs. (37) and (38), and 4>y can be expressed as the result of a Taylor 
serie s 


and 


8f 




X 


£ x ( e PX* e PY 


) + 


X 


3e 


X 


af 


X 


X 9e. 


X 


6y + 8(6 ) 


X 


(58) 


6V f Y e Py) 


5 

3e. 


at 


'X 3e. 


X 


6 y + 8(6 ) 


X 


(59) 


As previously stated, closed-form representations for f„ and f are not 

A Y 

known but second- order approximations are 


-e 


4> 


x 


X 


X 


2L e l n 


+ 1 + 


(60) 


l 2 / 2l e l n 

E, 
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(61) 


4> 


Y 




3 e X 


e 


Y 


with L- T the null length of the actuator and L. the distance between engine 
N e- 

centerline (thrust axis) and the gimbal actuator attach points. Previous analy- 
sis has shown that approximating engine rotation by gimbal angle is consistent 
with the other approximations. Furthermore, the Taylor series will be trun- 
cated, and the expressions for gimbal motion become 


O' 


G x 


84 > 


X 


8e 


X 


6„ + 


84> 


X 


X 9e. 


X 


X 


a. 


84 ). 


8e 


X 


84 ). 


6, _ + 


X 8e. 


X 


Taking the partial derivatives in Eqs, (62) and (63) yields 


“G„ ■ L„ [ (a l ' 11 6 x + a 2 6 y] 

X E 

and 

a G y = L^[ a 3 6 X + (a 4 “ 1} 6 y] 


(62) 


(63) 


(64) 


(65) 
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where 



Substituting the values 

L<£ = 6.43 cm (Z.53 in.) 

= 8.Z6 cm (3.Z5 in.) 

'L'N - 17. Z cm (6.763 in.) 

will yield the values of and a Gy as f unct i° ns of the pre-aim. For 

Mission B, the presently predicted values are listed in Table 5. Substitution 
of these values into Eqs. (66) through (69) will yield the coefficients for the 
gimbal motion Eqs. (64) and (65) and are listed in Table 5, 

Another source of error is that the gimbal reference axes will not be 
exactly aligned with the IRU sensed axes. Let the misalignment rotation be 
V radians so that 


16 


JPL Technical Memorandum 33-703 




cos 


sin 


V 

y 


— 


— — - 

sin y 


e x 

cos y 


e Y 



— * 


and via small angle approximations 


I 

CO 

X 

<D 

1 


1 

-y 


e x 

CD 

X 

CO 


V 

l 


®Y 

L 


, 




_ — 


( 70 ) 


(71) 


The 3 q- estimate for y is taken as 5 mrad. 

E« Steady-State Thrust Pointing Accuracy 

In Section IV- B it was shown that to satisfy the path guidance criterion, 
it is necessary to select the forward gains 


K 


X 


= K, 



with being the lever arm from the actuator clevis to engine center of 

rotation. Estimates of the 3cr errors in the various gains which make up 
and I< Y were developed in Section IV-C and are listed in Table 4. Applying 
these errors to Eqs. (45) and (46) and computing the RSS total gives a 3 cr value 
of 0.24 for K v and K v . To determine the total steady-state error, first let 

K x = L E (l + ei ) (72) 


and 


K 


Y 


L e (1 + e 2 > 


(73) 


so that combining with Eqs. (64), (65), and (71) gives 
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(74) 

Expanding Eq. (74) and neglecting second-order terms and higher yields 



and using the identity matrix [i] becomes 


(75) 



where the elements of K are 


K n ' 

6 1 ' a i 

(77) 

K 12 = 

- (a 2 + Y) 

(78) 

K 21 = 

- (a 3 - y) 

(79) 

K 22 = 

e 2 “ a 4 

(80) 
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Rewriting Eq. (76) as 


[— — 1 


[— —1 

-1 

1- — j 

e x 


I + K 


O' 

G x 

6 Y 




Gy 

_ - 


— 


— — 


yields again to first order 



(82) 


The vector form of the steady- state pointing error Eq. (5) is 
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Substituting Eq. (82) into (83) gives 
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r - — 

T 

- — 


— — 

*x 
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a r 

G x 

+ 

» A 

A x 
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— f 
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(83) 


(84) 


Equation (2) can now be used to produce 
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( 85 ) 
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* CM X 
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I - K 


a A 

A x 

<P Y 
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Ay 
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— — 


Equation (81) is now used to estimate the steady- state pointing error. Sepa- 
rating CM and angular errors and expanding yields for CM errors 


<?> CM X K 11 °CM X + K 2 1 a CM y 
^CM y = K 12 Q, CM x + K 22 °CM y 


and for angular errors 


”a x = (1 - K H )ff A x - K 21 “a y 


<P A = - K ,7 a A +(1- K,,) o 


12 A 


22' “A 


Y 


( 86 ) 

(87) 

( 88 ) 

(89) 


A CM offset error of O'-,, 

CM 


and an angular error of 


a 


divide according to 


sin (^) 

CM CM 

A 

CM y = a CM COS ^ 
Q A = a A sin ^ r ) 


a 


A 


Y 


<* A cos (r) 


where ft and r are assumed independent random variables uniformly dis- 
tributed 0 to 2tt. Recalling that e^ and e^ in Eqs. (77) through (80) are 
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gaussian random variables with a 3tr value of 0, 24, and that the total pointing 
error is 



(90) 


the steady- state CM offset and angular errors can now be computed. This 
was done using a Monte Carlo technique on the digital computer. The results 
are given in the form 



( M C * 3 " t c ' 1 “CM 


(91) 


f A ) = (M A ±3tr A J “A 
'ss 


(92) 


where 

M„, cr 4 mean and standard deviation of the steady- state pointing 
error due to a unit CM offset 


and 

M , cr . = mean and standard deviation of the steady-state pointing 

-A. A 

error due to a unit angular error 

The mean and 3cr values in Eqs, (91) and (92) are given in Table 6, 


V. TOTAL TV C SYSTEM POINTING ERROR 

A. Combining Transient and Steady-State Errors 

The AV v (t) functions necessary to compute the transient pointing error 
were derived in Section III-D and are given by Eqs. (31) and (32) and Table 2. 
The form of the TVC system transient pointing error is given by Eq. (9). 

For a unit CM offset error, let 


p c (t) 4 


AV xc< t) 

AV(t) 


(93) 
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and for a unit angular error 


P A (t) - 


AV (t) 
XA U/ 

AV(t) 


(94) 


Equations (93) and (04) can be solved in the digital computer using the complex 
arithmetic capability of Fortran. This was done, and the resulting time func- 
tions are depicted in Figs. 5 through 16. Figures 5 and 6 give P^t) f° r the 
first midcourse maneuver. They are identical except that Fig. 5 is an expan- 
sion of the first 25 seconds of Fig. 6. Similarly, Figs, 7 and 8 show P^(t) 
for the first midcourse maneuver. For the first orbit trim. Figs, 9 and 10 
give P (t) and Figs. 11 and 12 give P (t), and for the second orbit trim 
following lander separation, Figs. 13 and 14 give P^(t) and Figs, 15 and 16 

give P A^ a 

Since ideal path guidance is implicit, P^(t) approaches unity asymptoti- 
cally and P^Ct) approaches zero (see Section II). To include the steady-state 
errors computed in Section IV, it is only necessary to combine Eqs, (91) and 
(92) with (93) and (94) to give 


Y 


CM 


(t) = [P c (t) + M c ± 3<r c ]a CM 

(95) 

(t) = [P A (t)(M A *3«r A )]a A 

(96) 


The transient plus steady-state error for each flight condition can now be 
computed using Eqs, (95) and (96) in conjunction with Figs. 5 through 16 and 
the statistical data given in Table 6. 

B. Verification 

A complete nonlinear simulation of the TVC system has been written in 
CSSL^ III for the Univac 1108 computer and includes the following: 

(1) Six degrees -of-freedom 

(2) Simulation of actual TVC electronics 

^ Continuous Systems Simulation Language. 
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(3) Spacecraft structural flexibility using hybrid coordinates 

(4} Propellant slosh dynamics modeled using a mechanical pendulum 
analogy 

For direct comparison with Figs. 5, 9, and 13 f typical runs with an initial 
1.0-degree CM offset for the first midcourse maneuver, first orbit trim, and 
post-VLCA^ release are shown in Figs. 17, 18, and 19. In each run, the 
structural and propellant slosh damping are set at their 3 & lower limits. 
Comparison of these figures shows nearly exact agreement with the nonlinear 
simulation and thus verifies the accuracy of the foregoing analysis. The 
somewhat oscillatory behavior in the simulation results is due to the dis- 
turbances introduced by propellant slosh dynamics. 


VI. SUMMARY 

Given estimates of the total CM offset in degrees and the total angular 
error in degrees, the AV pointing error due to the thrust vector control sys- 
tem can be estimated by using Eqs, (95) and (96) and the statistical data in 
Table 6. For the first midcourse maneuver the results are: 


CM<‘> " [P C U) + 0- 071 ±0„13]* C m 

(97) 

V A (t) = [ P A (t) (1-00 =*= 0-21)] <* A 

(98) 


where P c (t) is given in Fig. 5 or 6 and P A (t) in Fig. 7 or 8. The (±) terms 
are the possible 3cr variations. For the first orbit trim, 


CM (t) = [P c (t)+ 0.068 ±0.13] a CM 

(99) 

V A (t) = [p A (t) (1.00 ± 0.21)] a A 

(100) 


where P^t) is given in Fig. 9 or 10 and F^t) in Fig. 11 or 12. 


2 


Viking launch capsule. 
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Finally, for the second orbit trim. 


Y CM (t) = 

+ 0, 082 ± 0. id 

(101) 

Y A (t) = 

[p A (t) (1.04 ± 0.22)] a A 

(102) 

where P^t) is given in Fig. 

13 or 14 and P.(t) in Fig. 15 or 16. 
A 

Alter- 


nately, P^(t) and P^(t) in the above equations may be solved on the com- 
puter in the manner described in Section V-A. The above equations may be 
used for Mission A or Mission B and for other maneuvers within each mis- 
sion phase, since the TVC system electrical parameters are identical for the 
two spacecraft, and differences in mass properties are small. 
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Table 1. TVC system parameters 


Parameters 

First midcourse 
maneuve r 

First orbit 
trim 

Second orbit 
trim 

Engine thrust T, 
N (lbf) 

1334. 47 (300. 0) 

1334. 47 (300. 0) 

1334.47 (300. 0) 

Spacecraft mass M, 
kg (slugs) 

3437. 59 (235. 55) 

2332. 54 (159. 83) 

1174. 96 (80. 51) 

TVC forward gain 

k g k fc k a 

3. 333 

3. 333 

1. 105 

Lever arm L, 
m (ft) 

1. 927 (6o 323) 

2. 38 (7. 82) 

1.37 (4. 506) 

Path guidance time 
constant T^, s 

4. 86 

4. 86 

00 

« 

Rate to position 
time constant T^, s 

2. 0 

2. 0 

2. 0 

Gimbal servo lag 
T A- s 

0. 15 

0. 15 

0. 15 

Spacecraft moment of 
inertia J, kg-m 2 
(slug-ft 2 ) 

6775. 0 (4997) 

5144. 0 (3794) 

1626.0 (1199) 

Path guidance loop 
gain KpQi V/V 

4. 328 

4. 328 

2. 101 

Torque constant K^, 
N-m (ft-lb/rad) 

2571. 85 (1896. 90) 

3180. 75 (2346. 00) 

1832. 79 (1351. 8) 


fo 

01 







Constant 

First 

midcourse maneuver 

First 
orbit trim 

Second 

orbit trim j 

K i 

1. 748 

2. 912 

8= 688 

K 2 

3. 102 

0. 303 

0. 731 

a 

0. 175 

0. 184 

0. 188 

b 

3. 298 

0 G 878 

1. 108 

z 

2. 425 - jl. 855 

1. 608 + jl. 050 

4. 710 + j8. 586 

w 

1. 254 + jl. 200 

2. 460 + j3 0 387 

2. 572 + jl. 269 

z 

2. 425 + jl. 855 

1. 608 - jl. 050 

4.710 - j8. 586 

w 

1. 254 - jl. 200 

2. 460 - j3„ 387 

2. 572 - jl. 269 












Table 3. Select resistors for setting gains 


Resistor 

Pitch 

Yaw 

R i 


909. 0 

R_ 


2. 37 

2 



R 3 

866, 0 


R 4 

2, 37 

-- 

R 6 

CV 

ro 

i— ■ 

« 

o 

210. 0 

R 8 

267. 0 

255. 0 

R 9 

no. o 

107. 0 

R i6 

806. 0 

806. 0 

R 1 7 

806, 0 

806. 0 

R 1 9 

10. 0 

10. 0 

R 

10. 0 

10. 0 
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1. All values in kft. 


2. Complete list in Ref. 3. 
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Table 4, Gains and tolerances 


Gain 

Value 

Tolerance 

k gy 

K FG l 

170. 5 V/rad 
1. 041 V/V 

±5. 5% 
±3. 0% 

Kfc h 

3. 141 V/V 

±3 0 58% 


2. 101 V/V 

±5. 0% 

KPG H 

4. 328 V/V 

±6. 27% 

k dv 

0. 068 mA/V 

±2. 0% 

k fb 

10. 97 cm/V 
(4. 317 in. /V) 

±2 0. 8% 
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Table 5. Pre-aim values and coefficients for gimbal motion equations 


Maneuver 
(Mission B) 

e PX, m 
(in. ) 

e PY, m 
(in. ) 

First 

midcourse 

0. 00215 
(0. 00847) 

-0. 00122 
(-0. 04802) 

First 
orbit trim 

0. 000287 
(0. 01130) 

-0. 000287 
(-0. 01130) 

Second 
orbit trim 

0. 00309 
(0. 12146) 

-0. 00359 
(-0. 01412) 


0. 00125 -0. 0188 


0. 018 


-0. 00553 


- 0 . 0211 


0.00167 -0.00443 -0.00131 


0. 0404 


0. 0596 




















'K G a G s+1) 


ig. 1. Engine and gimbal 


FORWARD 

COMPENSATOR 


GIMBAL 

SERVO 


K . 
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Fig. 2. Single-axis linearized TVC system block diagram 


COMPENSATOR AND 
PATH GUIDANCE 



Fig 3. TVC system signal flow path 
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TIME, s 

Fig. 5.. Pointing error due to CM offset, 
tirst midcourse maneuver, 0 to 25 s 



TIME, s 

Fig. 6 .. Pointing error due to CM offset, 
first midcourse maneuver, 0 to 160 s 
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Fig. 11. Pointing error due to angular 
error, first orbit trim, 0 to 2 5 s 


Fig. 12. Pointing error due to angular 
error, first orbit trim, 0 to 160 s 
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Fig. 1 

e rror, 


5 . Pointing errors due to angular 
post VLCA release, 0 to 25 s 


Fig, 16. Pointing errors due to angular 
error, post V LCA release, 0 to 160 s 



Fig, 17, AV pointing error, first mid- 
course maneuver, Mission B, 1-deg CM 
offset, structural Model VI-F 
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Fig. 18. AV pointing error, first orbit 
trim, Mission B, 1-deg CM offset, 
post- VLCA release, structural 
Model VI-F1B 



Fig. 19. AV pointing error, 
post-VLCA release 
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